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We derive the exact ground-state energy of the one-dimensional Ising model in random 
fields taking values h, and —h with general probabilities. The random-field Ising model 
on a ladder is also analyzed by showing its equivalence to the random-field Ising chain with 
field values ±/i and for h < J. The zero-temperature transger matrix is used to obtain the 
results. 
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1. Introduction 

The one-dimensional random-field Ising model is one of the simplest examples of systems 
with quenched disorder. In the long history of its intensive studies, a number of results have 
been published which give exact solutions in the ground state. Derrida et al 1 first derived the 
exact expression of the ground-state energy by a recursive calculation of the zero-temperature 
transfer matrix. They also obtained the exact ground-state energies of the one-dimensional 
random bond Ising model in a uniform field, the random bond Ising ladder and the 3x3 square 
lattice system. Farhi and Gutmann 2 obtained the exact spin-spin correlation function of the 
one-dimensional random-field Ising model. Their approach is based on the ground-state spin 
configuration associated with the general rule that dictates which regions of the random field 
configurations necessarily determine the spin configuration. This method is related to another 
method which calculates exactly the energy and the entropy of the random-bond Ising model 
at zero temperature. 3 The application of the method proposed by Derrida et al to error- 
correcting codes was performed by Dress et al. 4 They investigated a simple error-correcting 
code through the zero-temperature properties of a related spin model. Further generalization 
of the method proposed by Dress et al was carried out by Kadowaki et al. 5 They calculated 
exactly the ground-state energies of the ±J random bond model and the site-random model 
on strips of various widths. Studies of the effects of continuous random fields in the ground 
state related to hysteresis are also found in Refs. [6] [7] [8]. 

In the present paper, we derive the exact ground-state energy of the one-dimensional Ising 
model in random fields consisting three possible values ±/i and using a generalization of the 
methods in Refs. [4] and [5]. The solution is piecewise linear as a function of the random- field 
strength, a property shared with the case of random fields ±h. We also find that the one- 
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dimensional Ising model in random fields ±h and is equivalent to the Ising spin ladder in 
random fields ±h as long as h is smaller than the exchange interaction J. This observation 
leads to the exact expression of the ground-state energy of the ladder under random field, 
which represents the first exact solution for the random-field Ising ladder. 

The present paper is composed of seven sections. In §2, we introduce our model and explain 
the formulation of our method based on the zero-temperature transfer matrix. An important 
feature of the ground state of the one-dimensional Ising model in random fields is explained 
in §3. The explicit evaluation of the probability distribution of a relevant quantity and the 
exact ground-state energy are shown in §4 and §5. A useful property for diagonalization of 
the transition matrix is also given. In §6, we consider the Ising ladder in ±h random fields 
and show that this model is equivalent to the one-dimensional Ising model with three possible 
random fields ±h and if h < J. The exact ground-state energy of the ladder model is derived 
using the result for the one-dimensional case. Finally in §7, results are discussed. 

2. Ising spin chain in the general random field 



2.1 Model 



Our problem is defined by the Hamiltonian 

N N 

h = —j s n -is n — y~] h n s n , (i) 

n=l n=l 

where J is the ferromagnetic coupling constant, S n (= ±1) is the Ising spin variable and N + l 
is the number of spins. We use an open boundary condition. The random field at site n consists 
of three possible values with probabilities 

p for h n = +h 

P{K) = { q for hn = , p + q + r=l. (2) 

r for h n = —h 

Almost all of the previous studies concentrated on the case q = 0. The partition function of 
this model with spin Sn on the right edge fixed is 

N 

Zn(Sn) =EE E II exp{/3(J5 n _ 1 5„ + h n S n )}, (3) 

So Si S N - in =l 

where (3 is the inverse temperature. We can construct the partition function Zn+i{Sn+i) 
recursively as 

Z N +i{S N +i) = ^2 Z N (S N ) exp((3(JS N SN+i + h N+ iS N +i))- (4) 

Sjv 
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To simplify this expression, we rewrite (4) by using a 2 x 2 matrix. Let us introduce the 
notation = Zn(Sn = ±1)- The recursion relation (4) can then be written as 



7+ 
%N+1 



Tn+1 



where Tjv is the transfer matrix 
with z = e@. 



-J—hjy 



7+ 

N 



J N 



yj—hjy 



(5) 



(6) 



2.2 Recursion relations in the zero -temperature limit 

In the limit of very low temperatures, the partition function Zn(Sn) can be expressed in 
terms of the ground-state energy as 

Z N (S N ) ~ exp((3x N (S N )) = z x ^ Sn \ (7) 

where xn(Sn) is the ground-state energy with the edge spin state Sn- If we use the expressions 
xjv(+1) = xn and xn(— 1) = xjy + 2ajy, we have in the zero temperature limit 

(8) 



(9) 
(10) 







z Xn 


. Z n . 







Substituting eqn. (8) into eqn. (5), we obtain the following recursion relations, 

x n+ i = x n + J + h n+ i + 2 max(0, a n - J) 
= f( a n) - h n+ i, 
where f(a n ) is a piecewise linear function defined as 



/(Or. 



+J for a n > J 



for 



-J < a„, < J 



(11) 



—J for a n < —J 

As shown below, from the recursion relation (9), the average of the ground-state energy over 
random fields reduces to the average of max(0, a n — J) over the distribution of the energy 
difference a n . We note that recursion relations (9) and (10) always hold for any types of 
random fields, not necessarily of eqn. (2). 

2.3 Ground- state energy 

Using the recursion relation (9), we can express the ground-state energy per spin as 



lim — 



1 

1 



logZ N 



= — lim — {xjv + max(0, 2a7v)} 

N^oo N 
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1 



' JV-l 

lim -f- <! y~] (x n+ i - x n ) + x + max(0, 2a N ) 
N— >oo iv z — * 
I n=0 



= -J — (p — - 2 ((max(0, a — J)}) , (12) 

where a is the energy difference in the limit N — > oo and the double brackets ((•••)) mean 
the average over the probability distribution of the energy difference P{a). To obtain the last 
line of eqn. (12), we have invoked the self-averaging property of the energy. 

According to eqn. (12), we have to derive the probability distribution of the energy dif- 
ference P(a). The derivation of P{a) is given in §4. Before going into the detailed evaluation 
of P(a), we explain another important feature of the ground state of the present model. 

3. Spin configurations in the ground state 

In this section, we consider the spin configuration in the ground state by evaluating the 
energy for specific configuration of random fields. From this argument, we derive an important 
property of the Ising model in random fields. 

We first consider the case q = {h n = ±/i). If the random-field strength h is larger than 
2J, we can easily obtain the ground-state spin configuration because all spins are forced to 
be parallel to the random field. However, for h < 2J, the ground state has non-trivial spin 
configurations as discussed now. 

First let us consider the following random-field configuration which consists of the field 
+h at all but a single site, 

••• + + + + + - + + + + + •••, (13) 

where the symbol + represents the field +h and — refers to the field —h. All spins are parallel 
to the random field when h is larger than 2 J. However, if h < 2 J, all spins are up (S n = 1). 
To see this, we evaluate the energy difference between the two spin configurations, 

E 1 = +2 J -h + E + (14) 

E 2 = -2 J + h + E+ (15) 

AE = E 2 -E 1 = -4J + 2h, (16) 

where E\ is the energy for all spins parallel to the field, E 2 represents the energy of all-up 
spin state and E + means the energy resulting from all sites except for the location of field 
—h. According to eqn. (16), for h > 2J, the ground-state energy is Ei, and it is E 2 otherwise. 
We note that the ground-state energy is doubly degenerate at h = 2J. 

Next we consider the following configuration of the random field which consists of the field 
+h at all sites except for two adjacent sites 

■■■ + + + + + -- + + + + + ■■■. (17) 
Similar arguments as above lead to the conclusion that all spins point parallel to the random 
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field if h > 3, and all spins are up if h < 3. The ground-state energy is degenerate at h = 3. 

We can apply these arguments to the general case with k pieces of successive random fields 
taking the value —h, 

••• + + + + + 1 • + + + + ■••, (18) 

k 

for which 

E x = +23 -kh + E + (19) 

E 2 = -23 + kh + E + (20) 

AE = E 2 -E 1 = -A3 + 2kh, (21) 

where the definitions of E\, E 2 and E + remain the same as in the case of eqn. (13). The ground- 
state energy changes from E\ to E 2 at h = 23 /k and these two energies are degenerate at 
h = 23 1 k. Consequently, the ground state of the one-dimensional Ising model in the random 
field changes at the random-field strength h = 23 /k, where k is an integer. 

To understand it, we consider the following configuration of the random field, 

... + + +^+ + + ... + + +^+ + + --- + + +^-^+ + + -- - , (22) 
ABC 
where A, B and C are the clusters of fields pointing to the — direction with length one, two 
and three. There are several clusters of fields taking the value — h in the sea of fields taking 
the value +h. In the case of h > 23, all spins are parallel to random fields. When h decreases 
and satisfies h < 23, the spin at the cluster with length one (that is, the part A in (22)) flips 
its direction from down to up. Next, spins at the cluster of length two (the part B in (22)) 
flip their direction from down to up when h decreases and satisfies h < 3. The ground-state 
proceeds with this change, and spins in the cluster with length k flip their direction from —h 
to +h when the random-field strength reaches h = 23 /k. The ground-state spin configurations 
are invariant in the range 

--^—-3 <h<y3. (23) 
k+l k K ' 

Within this range the ground-state energy is a linear function of the random-field strength h. 
The same discussion as above holds if we include the possibility of random-field value of 0. 

4. Probability distribution 



For explicit evaluation of the ground-state energy, we should find the probability distri- 
bution of the energy difference P(a) according to eqn. (10). In this section, we calculate this 
probability distribution P(a) by the recursion relation (10) following the idea of Refs. [4] and 
[5]. 
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Let us see how the a's are generated stochastically by the recursion relation (10). For 
example, we first consider the simplest case of the random field strength h > 2 J. We start the 
recursion relation from the initial condition oq = and obtain a± = —h for hi = h, a± = for 
hi = and a\ = h for hi = —h. The next value of a n is 



02 


= -J -h 


for 


01 


= -h, 


h 2 


= h 


02 


= -J 


for 


ai 


= -h 


h 2 


= 


02 


= -J + h 


for 


ai 


= -h, 


h 2 


= -h 


0,2 


= -h 


for 


a\ 


= 0, 


h 2 


= h 


0,2 


= 


for 


CLl 


= 0, 


h 2 


= 


02 


= +h 


for 


CLl 


= 0, 


h 2 


= -h 


02 


= J-h 


for 


a\ 


= h, 


h 2 


= h 


02 


= J 


for 


CLl 


= h, 


h 2 


= 


02 


= J + h 


for 


CLl 


= h, 


h 2 


= -h 



(24) 



No other values emerge for a n (n > 3) as can be verified by repeating this procedure. Conse- 
quently, we can regard the recursion relation (10) as a stochastic process with the transition 
matrix 

r P + Q P P 

•S n Q (25) 

R R Q + R 



where 





p 










p 










p 


p = 


q 








, Q = 


q 





, R = 





q 




r 










r 










r 



(26) 



The row and column of matrix (25) correspond to the following vector of the energy differences 

-J-h 

-J 
-J + h 

- (27) 
+h 
J-h 
J 

J + h 

For example, the (1, l)-component of matrix (25) means that the energy difference —J — 
h remains unchanged with probability p. The stationary distribution P(a) is given by the 
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eigenvector of this transition matrix with the eigenvalue 1, and the result is 

P(-J-h) 

P{-J) 
P(-J + h) 
P(-h) 

P(0) - . (28) 

P(+h) 
P(J - h) 

P{J) 
P(J + h) 

This result holds in the region h > 2 J, the case of k = in eqn. (23). The ground state is 
invariant in the region of h > 2 J. This means that the stationary state of the range h > 2 J 
is always described by the transition matrix So and the probability distribution P(a) remains 
invariant throughout the range h > 2 J. A similar property holds in other ranges specified by 
eqn. (23). 

Applying the same argument to the case of k ^ 0, we obtain the following 3(2k + 3) x 
3(2k + 3) transition matrix 





" p 2 " 
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(29) 



The row and column of the matrix Sk correspond to the following vector of the energy differ- 
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(30) 



-J-h 
-J 

-J + h 
-J+h-h 

-J + h 
-J + h + h 

-J + kh-h 

-J + kh 
-J + kh + h 
-h 


+h 
J — kh — h 

J -kh 
J -kh + h 



J-h-h 

J-h 
J -h + h 
J-h 
J 

J + h 

Evaluation of the eigenvector of the matrix Sk is not simple because this matrix includes 3x3 
matrices P, Q and R. This difficulty can be resolved by the decomposition of the eigenvector 
into a direct product. The details are explained in Appendix. We only show the result here. 
The eigenvector of the matrix Sk with eigenvalue 1 is given as 



/ V \ 



V = Vq 



(31) 



V r J 
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where vq is the eigenvector of the (2k + 3) x (2k + 3) matrix 
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(32) 



with eigenvalue 1. 



4-2 Probability distribution 



Diagonalizing the transition matrix (29), we obtain the probability of the energy difference 
P(a). For this purpose, the relation (31) is useful. We define the following eigenvector of the 
matrix (32), 

v T = (£>i, D 2 ,- - ,D k ,D k+1 ^,D^D^--- (33) 
Using this expression of the vector vq, we obtain the following results 



fc+l / \ l-l k+l 
( p J 



l-i 



n - l=m ' ' n_ - ™ c\A\ 

i=i Vj y i=i Vj y i=i i=i 

where integer m runs from 1 to k + 1 and is specified by the range of the random-field strength 
as in eqn. (23). The probability distribution of the energy difference P(a) corresponds to the 
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P(—j -\-h- h) 




pis I 


P(— J + h) 




a Do 


p(- J + h + h) 






P( — J + Kfl — 11) 






r( — J + Ktl) 




qJ^k+i 


r\ — J -\- Ktl -\- 11) 
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pD T 
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rD T 



This is the final expression of the probability distribution of the energy difference P(a). 
5. Exact energy 



Finally, we obtain the ground-state energy for 2J/(k + 1) < h < 2J/k using eqn. (35) as 

f A{pr) k +\p-rf \ 

EGS = \- 1 + A(k + l)A(k + 2)S J 

/ , 2(pr) k + 1 (k + l)(p-r) 2 2A{k + l)pr \ 

\ q+ A(k + l)A{k + 2) A(k + 2) r' W 

where A{k) is 

A(k) =p k - r k . (37) 

In Fig. 1, we plot the ground-state energy eqn. (36) as a function of the random-field strength 
h for several probabilities. 
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-3.5 1 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 

h 

Fig. 1. The ground-state energy of the random-field Ising chain with h n = ±h and 0. 

6. Ising spin ladder in bimodal random fields 

Next, we consider the Ising spin ladder in bimodal random fields ±h. For this model, 
we can prove its equivalence to the one-dimensional model discussed in the previous sections 
when h < J. The exact ground-state energy of the ladder model is calculated exactly using 
the result in eqn. (36). 

6. 1 Spin configuration for h < J 

It is possible to show that for the random-field strength h smaller than J, the Ising spin 
ladder in the bimodal random field is equivalent to the one-dimensional Ising model in random 
fields ±h and 0. 

Let us consider the specific random-field configuration which consists of the fields being 
— h in the first chain and +h in the second chain (Fig. 2). For this random- field configuration, 
there are two possible ground states. One consists of all spins being up (or down due to the 
symmetry), and the other is all spins being parallel to the field. We introduce the notation 
-^aiiup for the energy of the first ground state and -Efidds fo r the energy of the second ground 
state. The energies -E a ii up and -Egeids are 

Gallup = — 2(n — 1) J -nh + n(h - J) (38) 

^fields = -2(n-l)J-nh + n(J-h), (39) 

where n is the number of the rungs. If the random-field strength h is smaller than J, the 
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+h +h +h +h +h +h +h +h 

Fig. 2. The configuration of fields being — h in the first chain and +h in the second chain. 
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A B 

Fig. 3. The general random-field configuration. The spins in parts A and B should be parallel for 
h<J. 



ground-state energy becomes -E a u up , and the ground state is all spins being up or down with 
equal probability. Let us consider a general random- field configuration as shown in Fig. 3. All 
spins in the parts A and B should be up or down from the above discussion. The spins in the 
left-most rung should be up because the fields on this rung are +h, and the spins between the 
parts A and B are also up for the same reason as the left-most rung. Accordingly, the spins in 
the part A should be up. The spins in the right-most rung should be down for the fields on this 
rung being down, and there are two degenerate states between all spins in the part B being up 
and down. The same argument as above holds for any other random-field configurations. For 
any case, the spins in the same rung are parallel. Consequently, we conclude that for h < J, 
two spins on the same rung in the ladder model always point to the same direction. Therefore, 
by identifying spins on the same rung as a spin taking the values ±1, we obtain the effective 
coupling constant 2J and the three kinds of effective fields {+2/i, 0, — 2h}. This means that 
the ladder model in the bimodal random field is equivalent to the Ising chain in random fields 
±h and 0. 
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6.2 Exact energy 

From above argument, the exact ground-state energy of the ladder model can be obtained 
from the result in eqn. (36). The corresponding probabilities of the random field are 

p = x 2 for +h 

q = 2x{l-x) for , (40) 

r = (1 - x) 2 for —h 

where x means the probability of the random field being +h for the ladder model. We also 

replace the random field strength and the coupling constant of the one-dimensional model 

with J — > 2J and h — > 2h. We need an additional term J which is caused by spins on the 

same rung pointing to the same direction, and divide the result by 2 because the number of 

spins in ladder model is twice as many as the one-dimensional model. The final expression is 

f 3 4(pr f+\p - r) 2 } 
GS 1 2 A(k + V)A(k + 2) J 

2{pr) k+ \k+l)(p-r) 2 2A{k + l)pr \ 
q A(k + l)A(k + 2) A(& + 2) J ' 1 j 



where -A(fc) is 



7. Discussion 



A(k) =p k - r k . (42) 



In the present paper, we have derived the exact ground-state energy of the random-field 
Ising chain with field values ±h and using the zero-temperature transfer matrix method in 
Refs. [4] and [5]. We have also obtained the exact ground-state energy of the Ising ladder in 
bimodal random fields for h < J. As far as we know this is the first example in which an 
exact solution has been derived for ladder model in random fields. The present techniques is 
a step toward a solution of random- field Ising models on strips with wider widths, eventually 
reaching the two-dimensional model. 

We give a few comments on some aspects other than the ground-state energy. First, it 
should be possible to derive the exact solution for other physical quantities than the ground- 
state energy, such as the entropy and magnetization, using the methods in Refs. [1] and [9]. 
Investigations in this direction are under way. 

For the Ising ladder in bimodal random fields, the evaluation of the transition matrix is 
not straightforward for h > J since we cannot map the problem to the one-dimensional model. 
This difficulty is caused by the complicated structure of the recursion relation of the ladder 
model. In fact, we find four kinds of recursion relations to evaluate the probability distribution 
of the energy differences. It is not easy to construct the general form of the transition matrix 
from these recursion relations. This is nevertheless an important and challenging task because 
of the eventual target of the two-dimensional model as mentioned above. 
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Appendix: Decomposition of eigenvectors into a direct product 

The matrix (29) is not simple to diagonalize because it is composed of the 3x3 matrices 
P, Q and R. To diagonalize this 3(2fc + 3) x 2(2k + 3) matrix, we use the decomposition of 
the eigenvector into a direct product. We define the following vector w, 

( P \ 



w 



(Al) 



V r J 

This vector satisfies the following relations 

Pw = pw, Qw = qw. (A-2) 

Using eqn. (A-2), we can decompose the eigenvector v of the matrix Sk into a direct product 
of vectors vq and w as 

v = vq <g) w, (A-3) 

where i>o is the eigenvector of matrix §k in eqn. (32) with eigenvalue 1. To prove eqn. (A-3), 
we show that the vector vq (g) w satisfies the equation for the eigenvector of the matrix Sk 
with eigenvalue 1 as 

SkV = Sfc(i>o <8> w) = {SkV ) <g> w = v <g) w = v. (A-4) 

Thus eqn. (A-3) holds. Accordingly, it is useful for the evaluation of the probability distribution 
that we calculate the eigenvector of the matrix Sk with eigenvalue 1 instead of the eigenvector 
of the matrix Sk with eigenvalue 1. 
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